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Using first-principles calculations, we identify "magic-zero" optical wavelengths, 
the ground-state frequency-dependent polarizabilities of alkali-metal atoms vanish. Our approach 
uses high-precision, relativistic all-order methods in which all single, double, and partial triple 
excitations of the Dirac-Fock wave functions are included to all orders of perturbation theory. We 
discuss the use of magic-zero wavelengths for sympathetic cooling in two-species mixtures of alkalis 
with group-II and other elements of interest. Special cases in which these wavelengths coincide with 
strong resonance transitions in a target system are identified. 
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I. INTRODUCTION 

The realization of mixtures of trapped ultracold 
atomic gases [ll-Q has opened new paths towards the for- 
mation of ultracold diatomic molecules [6|4lC| , quantum- 
state control of chemical reactions [ll[, prospec ts for 
quantum computing with polar molecules |12h14J . tests 
of fundamental symmetries [15l4l7l | and studies of fun- 
damental aspects of correlated many-body systems [18j . 
and dilute quantum degenerate systems |19l - [23| . Co- 
trapped diamagnetic-paramagnetic mixtures have also 
made possible experimental realization of interspecies 
Feshbach resonances |24l42a |. two-species Bose-Einstein 
condensates and mixed Bose-Fermi and Fermi-Fermi de- 
generate gases [il,!!^^. 

In an optical lattice, atoms can be trapped in the inten- 
sity maxima or minima of the light field by the optical 
dipole force [30|. This force arises from the dispersive 
interaction of the induced atomic dipole moment with 
the intensity gradient of the light field, and is propor- 
tional to the ac polarizability of the atom. When its ac 
polarizability vanishes, as can happen at certain wave- 
lengths, an atom experiences no dipole force and thus 
is unaffected by the presence of an optical lattice. Our 
present work provides accurate predictions of the "magic- 
zero wavelengths" Azero which lead to zero Stark shifts for 
alkali-metal atoms. These wavelengths have also been 
designated as "tune-out wavelengths" by LeBlanc and 
Thywissen |3l| . 

We suggest some possible uses for such wavelengths, 
all of which take advantage of the fact (demonstrated 
below), that magic-zero wavelengths are highly depen- 
dent upon atomic species and state. For a given atomic 
species and state A, let La designate an optical lattice or 
trap made with light at one of the magic-zero wavelengths 
of A. We start with a model configuration consisting of 



the gas A embedded in La and confined by another trap, 
T. Some process is performed on the gas, after which T 
is turned off. Members of A will depart and La may 
confine whatever is left. For example, one might pho- 
toassociate some A atoms into dimers during the initial 
period, and thereby be left with a nearly pure population 
of dimers trapped in La at the end. LcBlanc and Thy- 
wissen [31] have pointed out the advantage of magic-zero 
wavelengths for traps containing two species. If another 
species, B, is added to the model configuration, it will 
ordinarily be affected by La, so B can be moved by shift- 
ing La, while A remains unaffected. Schemes of this type 
have been used for entropy transfer and controlled colli- 
sions between ^^Rb and ^^K 32 34] . For bichromatic op- 
tical lattice schemes, such as those discussed by Brickman 
Soderberg, et al. [33, ISg , it could be useful to incorpo- 
rate Lb into the model configuration, so as to be able to 
move A and B completely independently. In another ap- 
plication, a Sr lattice at a ^Fq magic-zero wavelength was 
sug gest ed for realization of quantum information process- 
ing HJ- 

In the next section, we briefly discuss the calculation 
of frequency-dependent polarizabilities of alkali-metal 
atoms. In section IIIIl we present the magic-zero wave- 
lengths for the alkalis from Li to Cs and discuss some of 
their applications. 



II. FREQUENCY-DEPENDENT 
POLARIZABILITIES 



The background to our approach to calculation of 
atomic polarizabilities is treated in a recent review arti- 
cle [33] ■ Here we summarize points salient to the present 
work. The frequency-dependent scalar polarizability. 



Oq (w) , of an alkali-metal atom in its ground state v may 
be separated into a contribution from the core electrons, 
acoro, a core modification due to the valence electron, 
a^c, and a contribution from the valence electron, a^iuj). 
Since core electrons have excitation energies in the far- 
ultraviolet region of the spectrum, the core polarizabil- 
ity dependends weakly on lo for the optical frequencies 
treated here. Therefore, we approximate the core polar- 
izability by its dc value as calculated in the random-phase 
approximation (RPA) |39| . an approach that has been 
quite successful in previous applications. The core po- 
larizability is corrected for Pauli blocking of core- valence 
excitations by introducing an extra term a.uc- For con- 
sistency, this is also calculated in RPA. Therefore, the 
ground state polarizability may be separated as 



ao(a;) = acorc + a^c + ctl{uj). 



(1) 



The valence contribution to the static ac polarizability is 
calculated using the sum-over-states approach [40|: 
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where {k \\D\\ v) is the reduced electric-dipole (El) matrix 
element. In this equation, w is assumed to be at least sev- 
eral linewidths oflf resonance with the corresponding tran- 
sitions. Unless stated otherwise, we use the conventional 
system of atomic units, a.u., in which e,mo, 47reo and 
the reduced Planck constant h have the numerical value 
1. Polarizability in a.u. has the dimension of volume, 
and its numerical values presented here are expressed in 
units of Uq, where uq ~ 0.052918 nm is the Bohr radius. 
The atomic units for a can be converted to SI units via 
a/h [Hz/(V/m)2]=2.48832xl0-*a [a.u.], where the con- 
version coefficient is AneoaQ/h and the Planck constant h 
is factored out. 

The calculation of the ground state frequency- 
dependent polarizabilities in alkali-metal atoms has been 
previously discussed in |4ll . |42| , and we give only brief 
summary of the approach. The sum over intermediate 
k states in Eq. ^ converges rapidly. Therefore, we 
separate the valence state polarizability into two parts, 
Q^main, Containing the contributions from the few lowest 
np states, and the remainder, ataii- We note that our 
calculations are carried out with the finite basis set con- 
structed using B-splines [43| making the sum finite. In 
the calculation of amain) we use the experimental val- 
ues compiled in Ref. [4J| along with their uncertainties 
for the first ns — np matrix elements, for example the 
4s — 4pj matrix elements in K. For all other terms, we 
use the relativistic all-order values [4J: E^ of the matrix 
elements and the experimental values of the energies [461 - 
l48j . In the relativistic all-order method, all single-double 
(SD) or single-double and partial valence triple (SDpT) 
excitations of the Dirac-Fock (DF) wave function are in- 
cluded to all orders of perturbation theory [40, H, [43] ■ 
We conduct additional semi-empirical scaling of our all- 
order values (SDsc) where we expect scaled values to be 



TABLE I: 5s — np contributions to the frequency-dependent 
polarizability of the ground state of Rb at Azero = 
423.0448 nm = 1/(23638.16 cm"^). Absolute values of 
electric-dipole matrix elements are expressed in a.u. (eao), 
and the corresponding energy differences are expressed in con- 
ventional wavenumber units (cm^^ ). 



Contribution 



|(5s||D||npi/2)| 



En 



Es 



ao 



5pi/2 
Gpi/2 
7pi/2 

8Pl/2 

npi/2 

5P3/2 
6p3/2 
7P3/2 
8p3/2 

np3/2 
acore 

Total aaiuj) 



tail 



tail 



4.231 


12579.0 


-41.130 


0.325 
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0.115 
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FIG. 1: The frequency-dependent polarizability of the Rb 
ground state. The first two zero-magic wavelengths are 
marked with arrows. 



more accurate or for more accurate evaluation of the un- 
certainties. Our scaling procedure and evaluation of the 
uncertainties of the all-order results have been recently 
discussed in Ref. [5Q] . Briefiy, the uncertainties of the all- 
order matrix elements are given by the spread of their SD, 
SDpT, SDsc and SDpTsc values. These are also used to 
calculate the uncertainties in the state-by-state contri- 
butions to the frequency-dependent polarizability. The 
tail contributions, ataii, are calculated in the DF approx- 
imation using complete basis set functions that are linear 
combinations of B-splines [5l| . In the cases treated here, 
the tail contribution is of the order of 1% of the net va- 



TABLE II: Magic-zero wavelengths Azcro for alkali-metal 
atoms from Li to Cs. The resonant wavelengths Arcs for rele- 
vant transitions are also listed. The wavelengths (in vacuum) 
are given in nm. 



Atom Resonance 


Arcs 


^zcro 


''Li 2s - 2pi/2 


670.992478 


670.987445(1) 


2s - 2p3/2 


670.977380 




'Li 2s — 2pi/2 


670.976658 


670.971626(1) 


2s - 2p3/2 


670.961561 




Li 




324.18(6) 


2s - 3pi/2 


323.3576 




2s - 3p3/2 


323.3566 


274.911(7) 


2s - 4pi/2 


274.2001 




Na 3s — 3pi/2 


589.7558 


589.5565(3) 


3s - 3p3/2 


589.1583 


331.905(3) 


3s — 4pi/2 


330.3929 


330.3723 


3s - 4p3/2 


330.3319 


285.5817(8) 


3s — 5pi/2 


285.3850 




K 4s - 4pi/2 


770.1083 


768.971(3) 


4s — 4p3/2 


766.7009 


405.98(4) 


4s — 5pi/2 


404.8356 


404.72(4) 


4s - 5p3/2 


404.5285 


344.933(1) 


4s — 6pi/2 


344.8363 




Rb 5s — 5pi/2 


794.9789 


790.034(7) 


5s — 5p3/2 


780.2415 


423.05(8) 


5s — 6pi/2 


421.6726 


421.08(3) 


5s - 6p3/2 


420.2989 


359.42(3) 


5s — 7pi/2 


359.2593 




Cs 6s — 6pi/2 


894.5929 


880.25(4) 


6s — 6p3/2 


852.3472 


460.22(2) 


6s - 7pi/2 


459.4459 


457.31(3) 


6s - 7p3/2 


455.6557 


389.029(4) 


6s — 8pi/2 


388.9714 




lence contribution a^. 







We define the magic-zero wavelength Azero as the wave- 
length where the ac polarizability of the ground state 
vanishes. In practice, we calculated ao('^) for a range of 
values in the vicinity of relevant resonances and identi- 
fied values of uj where the polarizability turned to zero 



with sufficient numerical accuracy. 

We illustrate the cancellation of all contributions to 5 s 
Rb polarizability at Azero=423.0448 nm in TablelH Since 
this wavelength is between 5s — 5^3/2 and 5s — 6pi/2 
resonances, the contributions of the 5pj and 6pj terms 
strongly dominate. However, the contribution from the 
core is significant (11% of the largest valence term). This 
table shows that A^ero is located where the valence con- 
tribution to the polarizability cancels the adjusted core 
contribution, a feature that is common to all the cases 
treated here. The zero crossing point is in the close vicin- 
ity of the 5s — 6pi/2 resonance owing to the relative size 
of the 5s — 5pj and 5s — 6pj reduced electric-dipole ma- 
trix elements given in the second column of Table HI The 
5s — 6p matrix elements are more than an order of magni- 
tude smaller than the 5s — 5p matrix elements. Since po- 
larizability contributions are proportional to the square 
of the matrix element, the denominators of the 6pj terms 
have to become very small to cancel out the 5p contribu- 
tions. 

This magic-zero wavelength is illustrated in Fig. [T] 
where we plot ground-state polarizability of Rb atom in 
a.u. in the vicinity of the 5s — 6pj resonances. Another 
zero crossing point shown in the figure is located between 
5s — 6pi/2 and 5s — 6^3/2 resonances, as expected. The 
next magic-zero wavelength will be located close to the 
5s — 7pi/2 resonance since the values of the matrix ele- 
ments continue to decrease with n. 



III. RESULTS AND APPLICATIONS 



In Table |TT1 we list the vacuum Azoro wavelengths for 
alkali-metal atoms from Li to Cs. For convenience of 
presentation, we also list the resonant wavelength Arcs in 
vacuum in the relevant range of wavelengths. We order 
the lists of the resonant wavelengths and Azcro to indi- 
cate the respective placements of Azoro and their distances 
from resonances. The resonant vacuum wavelength val- 
ues are obtained from energy levels from National Insti- 
tute of Standards and Technology (NIST) database [4^ 
with the exception of the 2s — 2pi/2 and 2s — 2^3/2 tran- 
sition wavelengths for ^Li and ''Li that are taken from 
recent measurements |52l |. 

Since alkali ground states have electric dipole transi- 
tions only to p states, their polarizabilities will cross zero 
only between two ns — n'p resonances. We set the wave- 
length of the ns — {n + 2)pi/2 resonance as a lower wave- 
length bound for our search. The fine structure of the 
(n -I- 2)p level is sufficiently small for all alkalis to make 
the zero point between ns—{n+2)pi/2 and ns—{n+2)p^/2 
relatively difficult to use in practice, so we do not list it. 
We omit the Azcro between 2s — 3pj resonances for the 
same reason. The wavelengths of the next zero-crossing 
near the ns — {n + S)pi/2 resonances are in the ultravi- 
olet, and not as readily accessible in most laboratories, 
so we have not calculated them. However, this would be 



TABLE III: Wavelength (in vacuum) of selected transitions in Mg, Ca, Zn, Cd, Sr, Ba, Hg, Yb, Dy, Ho, and Er in nm. 



Atom 



Transition 



'So-3s3p 'Pi 
^So-4s4p Vi 
^So-5s5p 'Pi 
'So - 6s6p 'Pi 
'So-4s4p 'Pi 
'So - 5s5p 'Pi 
'So - 6sQp 'Pi 
'So ~ 6sQp 'Pi 
6s^ 5/8 - 4/1" (^/g) 6s6p i'Pi) J = 9 



Wavelength 


Transition 


Wavelength 


285.3 


3s^ 'So - 3s3p ''Pi 


457.2 


422.8 


As^ 'So - isip '■^Pi 


657.5 


460.9 


5s^ 'So - 5s5p ^Pi 


689.5 


553.7 


6s^ 'So'GsGp'^Pi 


791.4 


213.9 


4s^ 'So - 4s4p ^Pi 


307.7 


228.9 


5s^ ^50 - 5s5p ^Pi 


326.2 


184.9 


6s^ 'So - 6s6p ^Pi 


253.7 


398.9 


6s^ 'So-6s6p'^Pi 


555.8 


421.3 






582.8 






598.5 






416.4 







Mg 

Ca 

Sr 

Ba 

Zn 

Cd 

Hg 

Yb 

Dy 

Er 

Ho 

Ho 



3s" 

4s^ 

5s^ 

6s^ 

4s^ 

5s^ 

6s^ 

6s^ 
4/10 

4/ 
4/ 
4f 



'h 



5/2 



4f''{^He) 6s6p{'Pi) J = 7 

~ 4/11(4/15/2) 6s6p(3Pi) J =17/2 

- 4/11(4/15/2) 6s6p(iPi) J =13/2 



a routine matter for future work. There are no Azcro at 
wavelengths greater than those of the primary ns — npi/2 
resonances. Within these constraints, we have found four 
Azoro for Na, K, Rb, and Cs and three Azcro for Li, as 
shown in Table ID 

The stated uncertainties in the Azcro values are taken 
to be the maximum difference between the central value 
and the crossing of the ao i 6ao with zero, where Sao is 
the uncertainty in the ground state polarizability value 
at that wavelength. The uncertainties in the values of 
polarizabilities are obtained by adding uncertainties in 
the individual polarizability contributions in quadrature. 

We find small but significant differences in the first 
magic wavelengths of ^Li and ^Li due to the isotope shift. 
These values refer to the centers of gravity of all hyper- 
fine states and do not take into account the hyperfine 
structure. Therefore, this Azcro and the corresponding 
2s — 2pj resonance wavelengths are listed separately. We 
verified that isotope shift of the 2s — 2p traditions in Li 
does not affect the next magic wavelength, 324.18(6)nm, 
so we use NIST data for the other transitions. We also 
investigated possible dependence of the first magic-zero 
wavelengths on the isotope shift for '^^K, '^'^K, ''"'^K and 
**5Rb, 8^Rb. The Dl (4s - 4pi/2) and D2 (4s - 4^3/3) 
line wavelengths for '^^K, ^°K, and ^^K have been mea- 
sured using a femtosecond laser frequency comb by Falke 
et al. [53] ■ We carry out three calculations of the first 
magic-zero wavelength using Dl, D2 wavelengths for the 
specific isotope in our calculations. The resulting value 
for "^^K and ^"^K, 768.971 (3)nm, is the same as the value 
quoted in Table III The ^^K value is 768.970(3)nm, with 
the difference being well below our quoted uncertainty. 
The calculations of the first magic wavelength in Rb us- 
ing Dl and D2 frequencies for ^^Rb and ®^Rb listed in 



54 I55| gave results identical to result from Table |TI1 
790.034(7)nm that was obtained using NIST data. We 
note that our values for the first magic wavelen gth s are 
in good agreement with LeBlanc and Thywissen 3l| cal- 
culations with the exception of their value for *°K. 

The first magic-zero wavelength in Rb has been mea- 
sured to be 789.85(l)nm in [32i|. Some discrepancy with 



our result is most likely due to approximate linear polar- 
ization of the beam in |32l | . The difference is compatible 
with a shift in the magic zero wavelength caused by few 
percent spurious a~ polarization component |56| . 

Below, we identify two main applications of magic- 
zero wavelengths. First, these wavelengths are advanta- 
geous for cooling of group-II and other more complicated 
atoms, by sympathetic cooling using an accompanying 
alkali atom. 

Recently, group II atoms have been the subject of var- 
ious experiments and proposals in atomic clock research 
and quantum information. EEC of ®^Sr has been re- 
ported recently by two groups ^53, [5^ . The element Yb 
has four boson and two fermion isotopes, all of which 
have been cooled into the microkelvin range. Several ex- 
citing new prospects for quantum information processing 
with the ground state nuclear spi n have recently been 
identified in group II elements [33] ■ Also, Sr or Yb are 
useful for polarized mixtures of fermions, or Eose-Fermi 
mixtures, where isotopic mixtures can be studied. 

More complex systems have become of interest in the 
development of frequency standards and quantum infor- 
mation processing schemes. For example, the rare earth 
holmium is a candidate for quantum information appli- 
cations [5^ due to its rich ground-state hyperfine struc- 
ture. Erbium has been a subject of recent experimental 
work [60, [6l| , stimulated by its possible use in a variety 
of applications, including narrow linewidth laser cooling 
and spectroscopy, unique collision studies, and degener- 
ate bosonic and fermionic gases with long-range magnetic 
dipole coupling. 

Some species, particularly fermions, are difficult to cool 
by themselves due to unfavorable ultracold coUisional dy- 
namics. In such cases, it may be possible to use sympa- 
thetic cooling in a mixture of the target species and one 
of the alkalis, where the alkali-metal atom is cooled di- 
rectly by standard techniques. This has recently been 
demonstrated in Yb:Rb mixtures [62[. Use of Azcro trap 
wavelengths could allow one to release alkali atoms af- 
ter the target atoms of the other species are sufficiently 
cold, in a hybrid trap configuration that combines op- 



tical and magnetic traps or bichromatic optical traps. 
If the final trap configuration utilizes a A^ero wavelength, 
strong trapping of the target atom is possible while alkali 
atoms will be released by turning off its separate trap- 
ping potential. Since placement of the resonances varies 
significantly among the alkali-metal atoms, a wide range 
of Azcro is available as shown in Table IIIII 

We list the resonant wavelengths for variety of atomic 
systems in Table IIIII For consistency with the other ta- 
bles, we list vacuum wavelengths obtained from the NIST 
energy levels database [48| . Comparing Tables |ll] and Hill 
yields many instances of resonant transitions that are 
very close to Azcro- Here are few of the very close cases: 
Mg 285.3 - Na 285.6, Sr 460.9 - Cs 460.2, Dy 421.3 - 
Rb 421.1, Ho 598.5 - Na 589.6. 

Magic-zero wavelength laser light may be also useful in 
three-species cooling schemes such as reported in Ref. [6j| 
by allowing easy release of one of the species from the 
trap. The work [63] demonstrated that the efficiency of 
sympathetic cooling of the ^Li gas by ®^Rb was increased 
by the presence of '^'^K through catalytic cooling. 

Measurements of the magic-zero wavelengths may be 
used as high-precision benchmark tests of theory and de- 
termination of the excited-state matrix elements that are 
difficult to measure by other methods. Matrix elements 
of ns — n'p transitions of alkali-metal atoms, where ns 
is the ground state, are difficult to calculate accurately 
owing to large correlation corrections and small values 
of the final numbers. Experimental measurements of the 
-^zero predicted in this work will serve an an excellent 



benchmark test of the all-order calculations. Moreover, 
it will be possible to combine these measurements with 
theoretical calculations to infer the values of these small 
matrix elements. Only one high-precision measurement 
of such matrix elements (6s — 7pj transitions in Cs) has 
been carried out to date [6J|. 



IV. CONCLUSION 

In summary, we calculate magic-zero wavelengths in 
alkali-metal atoms from Li to Cs and estimate their un- 
certainties. Applications of these magic wavelengths to 
sympathetic cooling of group-II and other more compli- 
cated atoms with alkalis are discussed. Special cases 
where these wavelengths coincide with strong resonance 
transition in group-II atoms, Yb, Dy, Ho, and Er are 
identified. Measurements of the magic-zero wavelengths 
for benchmark tests of theory and experiment are pro- 
posed. 
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